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We argue that in clean high-zt type II superconductors, the low temperature vortex core size 
(defined as the coherence length ^) in high fields should decrease with increasing applied field in 
qualitative agreement with experimental data. Calculations are done for the Fermi sphere and 
cylinder (with the field parallel to the cylinder axis). The results for clean materials at T = can 
be represented as £,{H) /£^(Hc2) = U{H / Hc2) with U being an universal function. 
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I. INTRODUCTION 

The coherence length ^ has first been introduced as a 
phenomenological length scale in the near-Tc Ginzburg- 
Landau (GL) description where it sets, among other 
things, the "vortex core size" pc- Since the core, in fact, 
does not have a sharp boundary, the size pc cannot be 
unambiguously defined and is commonly chosen "opera- 
tionally convenient" , i.e., in a way which varies from one 
experimental or theoretical situation to another. For ex- 
ample, Pc may be defined as the radius of a circle where 
the persistent current is maximum,^ or - within the Lon- 
don approach - as the distance at which the divergent 
London current density reaches the depairing value, or 
assuming the core being a normal cylinder and equating 
the core contribution to the vortex energy to {H^ /8tt)ttPc 
with H^/Stt being the condensation energy.^''^ Another 
example comes from the scanning tunneling work where 
"the vortex core radius is arbitrary defined by that dis- 
tance Pc from the vortex center, for which the tunneling 
current has decreased from Imax to 36% of Imax~ Imin"'^ 
Clearly, these procedures yield different values of pc, al- 
though all of them have the same order of magnitude. 

A lot of experimental effort has been invested recently 
in study of the vortex core size, see the review by Sonier 
and references therein.^ Notably, whatever the definition 
of Pc is adopted, the low temperature pc is shown to 
decrease with increasing field in a number of materials 
such as NbSea, VaSi, LuNizBzC, YBazCugOr-i, CeRu2 
physical characteristics of which have little to do with 
each other (except all of them have a large GL parameter 
K — Xl/0- The dependences Pc{H) for all tested mate- 
rials are qualitatively similar; for large fields pc ~ 1/VH. 
Properties of the quasiparticle spectrum inside and out- 
side the cores (where the excitations may form narrow 
conducting bands) are considered responsible for the H 
and T dependences of Pc- 

The point of this paper is to provide a general 
theoretical argument based on the weak-coupling BCS 
theory in addressing causes of the field dependence of 
Pc- This argument is omitted in the current discussion 
of the problem in the experimental community, although 
a few examples of numerical solutions of the micro- 
scopic equations of superconductivity show that the H 
dependence of pc follows from the theory under very 



general assumptions.^"^ Our approach can be applied 
to the problem of Pc{H) only in large fields of high-K 
materials; still, we obtain our main results analytically 
that has certain advantages as compared to however 
powerful numerical techniques and enables one to make 
experimentally verifiable predictions. 

We begin with the notion that within the microscopic 
theory, at arbitrary magnetic fields H and temperatures 
T, it is not clear what exact value one should assign to ^. 
The difficulty comes from the fact that ^ and pc are not 
among the basic input parameters of the theory; instead, 
they should be calculated, analytically a very difficult if 
at all possible task in general. There is, however, a region 
at the second order phase transition from superconduct- 
ing to normal state, the SN boundary, where the theory 
can be linearized and consequently ^ is well defined. 

The linearization has been performed in the seminal 
work of Helfand and Werthamer (HW) on the upper criti- 
cal field Hc2{T, t) with r being the mean scattering time 
on nonmagnetic impurities.^ They have shown that at 
Hc2 where the order parameter A goes to zero, it satis- 
fies for any T a linear equation. 



-^2(T) n^A = A^ 



(1) 



where 11 = V + 27riA/0o is the gauge invariant gra- 
dient, A is the vector potential, and 4>o is the flux 
quantum. Formally, the equation is equivalent to the 
Schrodinger equation for a charge in uniform magnetic 
field; the field Hc2 at which superconductivity first nu- 
cleates corresponds to the lowest eigenvalue of this equa- 
tion: Hc2 — 0o/27r^^. This field (and £,) is obtained by 
solving the basic self-consistency equation of the theory: 



2^ ^ ^['^^ (3 - S 



where the function S{T, ^, r) can be written as: 

2/3 f"° , ,2 , s£q 
S = 4- dse^" tan^i — ^ 
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Here, uj = 7rT(2n + n is an integer, v is the Fermi 

velocity, and ^ = ut is the mean-free path. The power 
series representation of S is obtained by formally expand- 
ing tan~^ and then integrating over s. The evaluation of 
S was performed for the isotropic Fermi surface, i.e., for 
a Fermi sphere. 

Thus, strictly speaking, the length ^ is defined only at 
the SN phase boundary Hc^iT), and the question remains 
whether or not the same definition of ^ is useful out of 
the immediate vicinity of H(.2(T). In fact, in a variety 
of situations (small samples, proximity systems) the SN 
transition may take place far from the bulk Hc2{T). To 
approach the problem of the phase boundary in these 
systems, one has to know ^{H,T) in a broad domain of 
the H-T plane away of the bulk Hc2{T). 

A method to evaluate £^{H,T) had been developed in 
Refs. 9 and 10. In principle, the method follows HW 
by utilizing the field uniformity and the A smallness at 
the SN transition wherever it occurs. Below, we out- 
line the method as applied to the three-dimensional (3D) 
isotropic case of a Fermi sphere. Then, we consider 2D 
isotropic materials, i.e., the Fermi cylinder. We find that 
^{H,T) so obtained decreases when H increases toward 
Hc2 away of the critical temperature Tc in clean super- 
conductors; the effect is suppressed by impurity scatter- 
ing and is absent in the dirty limit. We provide a closed 
form equation for ^{H) for zero-T clean case for both 
the Fermi sphere and Fermi cylinder, and show that the 
results can be represented as 



u 



H 

Hc2 



(6) 



with U being an universal function. 

We next argue that the same procedure can be applied 
to the mixed state in applied fields H < Hc2 near the 
vortex core centers in materials with large k = XlH- 
This is because near the centers the field (varying on the 
scale of the London penetration depth Al) can be taken 
as uniform and A(r) at the center. We find our 
results in qualitative agreement with the data available, 
uncertainties of experimental procedures of extracting ^ 
notwithstanding. 



II. FERMI SPHERE: iTc2(T) AND S{H,i,u,t) 

Here we reproduce major points of the ^ [H) derivation 
for the system near the second order SN transition with 
the help of the quasiclassical Eilenberger formalism. 
The main equations of the theory read: 



TVUf 

A 



:r7(F + 2TA) 



2ttT 



,7; f^/A 

w=0 ^ 



(G 
F 



2u;T)f, 



(7) 
(8) 



Here, v is the Fermi velocity; /(r, w, v) and g{r, co, v) are 

the Eilenberger Green's functions with averages over the 
Fermi surface denoted as F = (/) and G = (<?). 



In the normal phase / = and g = I. In a small 
vicinity of the SN transition, |/| <C 1, whereas g can 
still be set unity in linear approximation in / due to 
normalization g = {1 — ff^Y^^ (for the same reason we 
do not need here an equation for f^). Equation (7) can 
be linearized: 

e-Uf = F-f3f, (9) 
i = VT, F = F + 2TA/h, /3 = 1 + 2a;T. (10) 

The solution of Eq. (9) is written as 

f = {P + e-U)-^F = r dpe-P^f'+^-^^F, (11) 
Jo 



or for the Fermi surface average: 
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dpe-Pf^ie-P^-^'F). 



(12) 



We now assume that A,F, and F satisfy Eq. (1); then, 
utilizing commutators of the operator II in an uniform 
field and the known properties of exponential operators,^ 
one can manipulate Eq. (12) to 



F{r,uj) = A(r) 



2tS 



(13) 



where 

oo 



m,j=0 



j!(2m + 2j + l) \ m\ J V/?^ 



2ttH 



(14) 



After substituting F of Eq. (13) in the self-consistency 
Eq. (8) and cancelling A(r), we obtain an implicit Eq. (2) 
for ^{H,T,£). It is easy to see that at J?c2 where = 
1/^^, the scries (14) reduces to the HW sum (4). 

The double sum (14) is in fact an asymptotic series and 
is diSicult to deal with when the goal is to solve Eq. (2) 
for ^ {H, T, £) . The situation simplifies greatly in the dirty 
limit: S is an expansion in powers of £. Keeping only the 
terms with m + j = 0, 1, one obtains 



1 - 



V?* < 1 , < 1 • (15) 



When substituted in the self-consistency Eq. (2), this 
yields de Gennes-Maki dirty limit result for ^(T, £) and 
Hc2{T)-^ Since the field does not enter S, the coher- 
ence length in the dirty limit is field independent. In other 
words, in the dirty limit, the coherence length at a given 
T determined at the upper critical field, is the same at 
this T for any H . 

Formally similar situation takes place near the critical 
temperature where the truncation (15) is justified by 
smallness of and ^"^ (not of £). We conclude that 
near the coherence length is field independent for any 
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t. We then expect the strongest H dependence of ^ to 
exist at low temperatures in clean materials. 

Given the complexity of series (14), it is desirable to 
have an integral representation for S better suited for 
analytic and numerical work. This had been done in 
Ref. 10 for the 3D case of a spherical Fermi surface. We 
refer the reader for details of this nontrivial procedure 
and provide here the result: 



S{u, a) — -^7? Re 







ds ^,. erfc5, (16) 



■2 ' 



2 V^^CV 



(17) 



erfcs = (2/y^) /^°^ c?z cxp(— z^) (cr here differs by the 
sign from that used in Rcf. 10). Note that ct = 1 at Hc2'; 
integration by parts in (16) gives the HW integral (3). 
One can check by formally expanding the integrand in 
powers of us^ and integrating over s, that the integral 
(16) can indeed be written as the series (14), see Ap- 
pendix B in Ref. 10. 



III. FERMI CYLINDER 

The calculations of the previous section cannot be done 
for an arbitrary Fermi surface. Still, for some simple 
shapes it is possible. The simplest of those is the Fermi 
cylinder with the field parallel to the cylinder axis. Em- 
ploying the same procedure outlined for the Fermi sphere, 
we arrive at 
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(-l)"2"(2m + 2n)! /mn 



n!(m!)^ 



(-) {l-a)m, (18) 



where u and a are defined in Eq. (17); we use the notation 
(1 - (j)m = (1 - o-)(2 - ct) • • • (to - The integral 

representation of this sum can be obtained in a manner 
similar to that described in Rcf. 10 for the 3D case: 



S{u, a) = —;=Re ds ^7- — e " 

see Appendix A. To verify this result we write: 

(1 - S^uf-^ _ /I - S^M^rr-l 1 



(19) 



1 + s'^uJ 1 + s'^u 

,2 



(1 + s'^uY 

a - 1 \ / _ _2us^\ 1_ 

\ m )\ l + us'^J 1 + 

m 

where fi = m + n. Substituting this to Eq. (19) and 
integrating over s one obtains 5 in the form (18). 

Having the quantity S{H, ^, T, £) for both 3D (Fermi 
sphere) and 2D (Fermi cylinder) , one can solve Eq. (2) 
for ^{H, T, £) that in general can be done numerically. As 
pointed out above, the most interesting is the situation 
at T = in clean materials; this case can be treated 
analytically. 



IV. CLEAN MATERIALS AT T = 

To deal with the divergence of ln(rc/r) in Eq. (2) we 
note than the sum over lu on the right-hand side (RHS) 
is actually extended to the Debye frequency LjJd- Then, 
we have for the finite sum 



u>0 



Hlo 2-kT 



In 



(21) 



where the neglected terms arc of the order T'^/h'^Luj^ and 
7 « 0.577 is the Euler constant. Hcucc the divergent InT 
in Eq. (2) drops off. Since in the clean limit /3 « 2ojt and 
2tS/{i3 — S) S/u) we obtain instead of Eq. (2) in the 
zero-T limit: 



In 



2hu!D 



2nTy^^ n ^S{u,a). (22) 
^of^ Jo ^ ' ' 



The integral at the RHS diverges logarithmically for 

ujd ^ oo, and so does the LHS; in other words, tj/) 
should drop off the result. This integral is evaluated in 
Appendix B for both 3D and 2D cases. 
We then obtain an implicit equation for ^: 



In 



hvq cos(7ri7) 
Aoa 4 



0, 



(23) 



where tp is the Digamma function. The only difference 
between 2D and 3D situations is in the number a: 



a2D = \/2 , a^D ^ e/V2. 
Setting cr = 1 in Eq. (23) one obtains: 

hv 



Hc2 = 



Cc2 = 



-7/2 



This yields for the 3D case 

/?c2(0) 



(24) 



(25) 



(26) 



the value obtained variationally by Gor'kov^'' and proven 
to be exact by HW; in HW reduced units it corre- 
sponds to ft,*(0) = H,2{^)/TcH'^^{Tc) « 0.72.8 p^j. ^-^e 
2D case, this gives /i*(0) w 0.59, the result obtained by 
Bulaevskii.^'* 

We now observe that material parameters enter 
Eq. (23) only in the first term under the log-sign. If one 

measures the length in units of ^c2 and uses the reduced 
field h = H/Hc2, Eq. (23) takes the form independent of 
material parameters: 



ln(2/ie^) + 
1 

"=2 



cos(7r(7) 



1 



h£,* 



?c2 



V(ct) =0, 
(27) 



4 



Hence, this equation defines an universal curve £,*{h) in- 
dependent of either material characteristics vp, Aq or the 
dimensionality. Given this curve and -ffc2(0), one can re- 
cover ^(i?) for a clean material at T = 0. 

The curve ^{H)/^{Hc2) = U{H/Hc2) is shown as a 
solid line in Fig. 1 for 0.15 < H/Hc2 < 1; the reason why 
the small fields domain is not shown is given in the next 
section. Also shown are results of numerical evaluation 
of ^ {H) for a few values of the impurity parameter A = 
hv/2TTTci. The numerical calculation is done with the 
help of the self-consistency Eq. (2) for arbitrary T and 
A; S{^, q^, T, A) is evaluated using an explicitly real form 
given in Appendix C. 

It is worth noting that effect of raising temperature on 
^{H) is qualitatively similar to that of the impurity scat- 
tering, see solid dots for t = T/Tc = 0.5: both suppress 
the field dependence of ^. However, at low temperatures 
for reasonably clean materials in a broad domain of high 
fields, ^{H) is well represented by the zero-T clean-limit 
curve; it is seen from the upper panel of Fig. 1 that for 
A = 0.25 this domain extends down to /i ~ 0.4. 

The lower panel of Fig. 1 shows the same results plot- 
ted against 1 / \/7i, the quantity proportional to the inter- 
vortex spacing. In this manner the data are often pre- 
sented to examine possible connection between the field 
dependence of the core size Pc{H) and other properties 
of the mixed state. ^ Our result shows that for materials 
on the clean side with A < 1, the slope d^* /d{h~^^^) for 
H Hc2 is universal. In fact, using Eq. (27) this slope 
at Hf.2 can be evaluated for the clean limit at T = : 

,,5i/2, =1-4-0-189. (28) 

We also observe that for real materials with A 7^ and 
T ^ 0, £,*{h) becomes flat as the field decreases with the 
impurity and temperature dependent plateaus. 
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FIG. 1: (Color online) The upper panel: the normalized co- 
herence length ^* = i{H)/^{Hc2) versus h = H/H^2- The 
solid line is calculated with the help of Eq. (27) for the clean 
limit at T = 0. The open simbols show ^*{X,t,h) for a few 
values of the scattering parameter A — hv/2TiTcl and reduced 
temperatures t = 0.1 shown in pairs {\,t) on the legend. The 
full symbols are for a clean material at an elevated tempera- 
ture: {\,t) = (0.1,0.5). The lower panel: ^* versus l/y/h. 



V. THE CORE SIZE 

The above discussion of £^{H) applies at the SN second 
order phase transition where the field is uniform and the 
order parameter A goes to zero. In fact, these conditions 
are met in vortex cores of high-K type-II superconductors 
in high fields. Indeed, in this case the field within the core 
of a size ^ is practically uniform since it varies on a much 
larger scale A^. Besides, when one approaches the vortex 
center, A ^ 0. To evaluate ^ in this situation, one can 
use the same formalism as at the SN phase boundary; 
in other words, the above procedure of evaluating £,{H) 
can be used to characterize the size pc- The core size so 
defined is, of course, "operationally convenient" from our 
point of view just as various definitions mentioned in the 
Introduction. Our approach, however, is advantageous 
because along with evaluation of the H dependence of 
the "core size" ^, we predict that 

• this dependence is weakened by scattering and dis- 
appears in the dirty limit, 



• the H dependence of ^ vanishes as T ^ T^, 

• ^ [H] is weakly affected by peculiarities of the Fermi 
surface, i.e., we expect qualitatively the same de- 
pendence for various materials, 

• in reduced variables, the dimensionless coherence 
length ^* = ^/^c2 should be nearly universal func- 
tion of the reduced field h = H/ Hc2 for clean ma- 
terials in high fields and low temperatures. 

Therefore, our "theoretically convenient" definition of 
the core size can be checked experimentally. 

One should put a note of caution on our claim of "uni- 
versality" as stated in the last two points. This feature 
expressed in Eq. (27), has been derived for two simple 
Fermi surfaces, a sphere and a cylinder. Nevertheless, 
since the Fermi surface shape always enters calculations 
of macroscopic parameters as ^ or Al via averaging over 
the whole surface, one does not expect the fine features 
of the surface to alter drastically our conclusion (except 
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in special circumstances, e.g., when the local density 
of states has sharp maxima at the surface). With this 
caveat we will use the term "universality" in further dis- 
cussion. 

There is another drawback to our approach. When 
applied to the SN phase boundary, say, of a proximity 
sandwich, the field H in the $(i?) dependence is the ex- 
ternally controlled uniform applied field. Defining the 
core size as £,{H), we imply that H is a the field value 
at the vortex center, Hq, which is nearly constant within 
the core provided k » 1. However, the problem is that 
there is no reliable and generally applicable estimate of 
Hq, except numerical results with a particular choice of 
parameters for 1ow-k and for high-T^ materials. An ex- 
ception is the case of isolated vortices (the applied field 
Ha 0) in the GL domain, where i/g ~ 2iJci-^^ In a 
more interesting situation of Ha ^ Hci , the vortex fields 
are strongly overlapped, and variations of the actual field 
within the vortex lattice are weak relative to the applied 
field; in fact, they are on the order of Hci <C Ha- There- 
fore, an error made by considering the field at the vortex 
axes as equal to the applied field is small. In other words, 
defining the vortex core size Pc{Ha) as £,(Ha) has a rea- 
sonable chance of success only in large fields Ha » Hd 
and improves as Ha — > i?c2- 

There is also a theoretical difficulty we encounter at- 
tempting to extend the analysis to /i ^ 0. In fact, the 
curve shown in Fig. 1 and generated by solving Eq. (27) 
shows oscillating behavior if extended to low fields be- 
yond h « O.IS.^'^ Our numerical work for finite A and T 
shows that these oscillations arc washed out quickly with 
increasing scattering and/or temperature and are hardly 
seen for A > 0.25, i.e., in still rather clean materials. 



VI. COMPARISON WITH DATA 

Relating the results obtained to information available 
on the vortex core size, we focus on the /xSR data re- 
viewed recently by Sonier."'^ This technique allows one to 
obtain the field distribution h(r) within the vortex lat- 
tice. Then one can calculate the current distribution and 
define the core radius pc as the distance from the vortex 
axis to the current maximum in the nearest neighbor di- 
rection. This definition of the core size is independent 
of a model one may choose to theoretically describe the 
distribution h{r), the point stressed in Ref. 18. We will 
consider here the data on so defined pc- 

The fxSR data on h{r) can also be analyzed with the 
help of the London model or its nonlocal version. For sim- 
plicity, we consider here the standard London isotropic 
result: 




where B is the magnetic induction and the sum is ex- 
tended over the reciprocal lattice G. 



It is relevant for this discussion that (a) the London 
model contains only one length scale, the penetration 
depth Xl, and (b) the model imphes the constant order 
parameter A and therefore breaks down at distances of 
the order ^. The latter comes about formally in Eq. (29) 
since the sum is divergent (this is readily seen as the 
logarithmic divergence of h when r — > 0). To mend 
this inherent shortcoming of the London model, various 
cutoffs arc commonly used, e.g., by introducing a fac- 
tor exp(— const G^^^) which excludes distances smaller 
than ^. Numerous efforts to fix the constant's value 
notwithstanding (see, e.g., Ref 19), in practice this con- 
stant is used quite liberally depending on the applica- 
tion in question. Other cutoffs basically suffer of similar 
uncertainties. Hence, the reliable results of the London 
model are only those that are insensitive to the cutoff 
chosen. Still, one can fit the data h{r) to a properly 
truncated sum (29), and extract the best-fit parameters 
Xl and ^ along with their H dependence. Interestingly 
enough, the so extracted £,{H) behaves as a function of 
H in nearly the same manner as Pc{H) extracted directly 
from the field distribution; it is found for a few materials 
that in high fields Pc^ ^ + C with a material dependent 
constant 

We now consider the data on Pc{H) for VaSi, NbSe2, 
YNi2B2C, and CeRu2 provided in Refs. 18,21- 23, respec- 
tively, and summed up in the review by Sonier.^ All 
the samples are high quality single crystals and have 
large GL parameters n: we assume them clean (the 
available scattering parameters are X{V3Si) « 0.13 and 
X{NbSe2) « 0.15). The reduced temperatures of the 
nSR experiments were low: « 0.22, 0.33, 0.19, and 0.3, 
respectively. For each material we have taken the Hc2 at 
a corresponding temperature, calculated ^c2, and normal- 
ized the experimental core radius to this value to obtain 
Pc = Pc/Cc2 • The results are plotted in the upper panel 
of Fig. 2 together with the theoretical versus reduced 
fields h = H/Hc2- For reasons explained above we took 
only the data points for h > 0.15. We expect the exper- 
imental p'*{h) and the theoretical £,*{h) to be shifted by 
a material, temperature and purity dependent constant: 
dh) ~ + Since the temperatures and 

impurity parameters in different experiments were differ- 
ent, we do not expect these shifts to be the same for the 
materials examined. In this situation, we have chosen 
the constants C so as to shift the data points as close 
as possible to our curve of CW- The result is shown in 
the middle panel of Fig. 2; the shifts needed are shown in 
the panel legend. Although the shifts vary, the data for 
different materials land nicely in a vicinity of our curve. 
This supports our guess of "universality" , a considerable 
ambiguity of the procedure notwithstanding. 

An interesting feature of the data and of the universal 
curve ^*(l/\//i) is seen in the lower panel of Fig. 2: the 
slope of this curve starting with the value (28) of « 0.2 
at /i = 1, increases to about 0.4 in the domain 0.25 < 
h < 0.5, and then drops back to about 0.2 near h » 0.15. 
In other words, the slope does not change much in each 
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FIG. 2: (Color online) The upper panel: the experimental 
core radius pc normalized on ^c2 = ^(t}o/2-KHc2{T) of each 
compound for materials indicated in the legend; T is the tem- 
perature of each experiment. The solid line is the theoretical 
with h = H/Hc2{T), the same curve as in Fig. 1. The 
middle panel: the same data shifted by amounts indicated in 
the legend. The lower panel: the same as the middle panel, 
but plotted versus l/\/7i. 



of these broad domains: 



dil/Vh) 



0.2 - 0.4, 



(30) 



or in common units: 



dil/VH) 



(0.2 




(31) 



Since the measured core size pc differs from our ^ by a 
constant shift, we may rewrite the last estimate as 



dpc 



d{l/VH) 



(115 - 230) AVkOe 



(32) 



in units employed in experiments. Given our sugges- 
tion of "universality" , we expect the high-field slope 
dpcl d(l/ ^/H) for all materials to be in this range. For 
the set of materials discussed here, this is the case: we 
roughly estimate the slope as 2 20 fo r VsSi, 190 for NbSe2, 
240 for YNizBaC, and 170 A%/fcOi for CeRua- The abil- 
ity of our approach to provide the slope values in a good 
agreement with the data indicates that the model catches 
correctly the physics of the field dependence of the core 
size. The slopes are relevant in particular, given an un- 
certain relationship between ^ we calculate and experi- 
mental Pc (uncertain shifts C in ^* {h) « p* (h) + C(A, T) 
shown in the middle panel of Fig. 2). 

Still, a number of questions remains to be addressed. 
Theoretically, it is not clear whether or not our clean limit 
results are compatible with the prediction of Kramer and 
Pesh that the core size of an isolated vortex defined as 
pi = A/ {dA/dr)r^o goes to zero as T ^ 0.^'* We note, 
however, that our results for pc are meaningful only in 
large fields and for large GL parameters k, whereas these 
authors have considered an isolated vortex in a mate- 
rial with K = 0.9 . The same can be said with respect 
to calculations of Ichioka et al. done for the d-wave 
symmetry who find a shrinking core size in decreasing 
temperatures.^^ 

Calculations of Miranovic et al. of the low temperature 
field dependence of the length pi = A„i/(9A/5r)r-+o in 
the mixed state (A^ is the order parameter maximum 
along the nearest neighbor direction) show that for clean 
materials with A < 1 the length pi{H) goes through a 
minimum and increases approaching Hc2- It also shows 
a much stronger field dependence on the dirty side (A > 
1) than our £,{H). A way out of this difficulty, in our 
opinion, is to conclude that pi{H) is not proportional 
either to our £,{H) or to existing data on the core size 
Pc{H) in clean materials for which the minimum in Pc{H) 
had not been recorded.^ 

We have to mention that our claim that the field depen- 
dence of ^ and Pc should be suppressed by impurity scat- 
tering and disappear altogether in the dirty limit contra- 
dicts calculations of Golubov and Hartmann done within 
the dirty limit Usadel formalism: they do find the H de- 
pendence of Pc for NbSe2.'* Their conclusion have been 
questioned in Ref. 21 on the grounds that the dirty limit 
approximation does not hold for NbSea crystals. More- 
over, Nohara et al. in Ref. 22 report that the H depen- 
dence of Pc extracted from the field dependence of the 
specific heat coefficient 7 and well pronounced in pure 
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NbSe2, in fact disappears after doping the crystal with 
Ta. The doping changes the impurity parameter from 
0.19 to 1.25 so that the observation is consistent with 
our conclusion. The data of this group on YNi2B2C and 
Y(Nio.8Pto. 2)620 are more convincing yet: the first crys- 
tal has the impurity parameter A w 0.4 (i.e., it is on the 
clean side) whereas the Pt-doped crystal is on the dirty 
side with A w 2.4; the H dependence of pc in the doped 
crystal is practically absent. Still, the controversy re- 
mains and could be resolved if the /xSR data were taken 
on a set of the same crystals with varying mean free path; 
candidates for such a study could be, e.g, the Co-doped 
LuNi2B2C crystals.26 
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APPENDIX A 

The sum (18) over m, n can be replaced with the sum 
over m from to ji = m + n and the sum over ji from 
to 00; the former can be written as a hypergeometric 
function: 

5 = ^(2m-i)!!(-|) 2Fi(-M,i-^;i;2). (Al) 

We now use the integral representation 



2Fi{a,b;c;z) = 



g— iTrfc^l — c 



r(c)r(i -I- 6 - c)r(i - 6) 



47r2 

dt{t- zf-^-H^-^{l-ty , (A2) 



where the contour circles the branch points at f = 
and t = z twice in opposite directions; the representa- 
tion holds everywhere except points where the F-factors 
diverge.^'' This yields: 

47r2 



-r(i - t7)r(t7) (b{t-2y-H 



J2{2i,-l)u\-^{l-t)Ydt. (A3) 



The sum here is transformed to an integral with the help 
of identity 



^(2M-l)!!(-a.)'' = -^Re ' 



1 -I- 2s^x 



(A4) 



which is proven by formally expanding 1/(1 -|- 2s^x) in 
powers of 2s'^x and integrating over s. 

Then the contour integral emerges of the form 



a—li—a 



J= (b{t-2f-H 



(It 



1 + s^u{l - t) 



(A5) 



which can be transformed back to the hypergeometric 
form after the substitution t = v{l + s^u)/ s^u = v z: 



47r2gj'r(i-<^) 

r(i - fT)r(a) 



2Fi(l,l-c7;l;2/^). (A6) 



Further, 2^1(1, 1-ct;1; 2/2) = (l-2/0)''-i and we obtain 
Eq. (19) of the main text. 



APPENDIX B 

One has to evaluate the integral at the RHS of Eq.(22); 
we start with the Fermi cylinder: 

/ — S = —=^% \ dse — r]((j,uj), (Bl) 

Jo Vt^ Jo Jo ^ 



with 



2„,^<T-l 



(1 - s^u) 
{1 + s'^uy 



u = 



4^ 



(B2) 



Substitution x = w^/wf, transforms the integral over uj 
to 



= !/■ 

2 70 



dx 



{x - vY 



{x + vY 

{-l)-{l + yY- 

' 4y((7-l)(2y)- ' 



Fi 1 - cr, 1 - cr; 2 - cr; 



2y 



where 



(B3) 



(B4) 



because large values of s are cut off by e^^ . Utiliz- 
ing the reflection formulas for the Digamma function,^^ 
V'(l — z) = tp^z) + 7rcot(7r^), the expression in square 
parentheses is rewritten as 



27r 



sin(7rcr) 



(B5) 



We further use the asymptotic formula 15.3.13 of Ref. 12 
for 2F1 since l/2y 1. Then, the real part of J assumes 
the form: 



Re J 



cos(7rcr) 



H^)-Hf) 



[ln2-|-7-F'^(a-)] -l-ln— ^ - Ins. 



(B6) 
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The integration over s in Eq. (Bl) is now straightfor- 
ward: the s independent part of ReJ enters the result 
being unchanged because dse~^ = 1. Fur- 

ther: (2/V^)/Q°°d.se-'''lns = -7/2 - ln2. Collecting 
all terms in the self-consistency equation (22) we obtain 
Eq. (23). As expected, the large parameter ivn cancels 
out from the final result. 

For the 3D situation, we have to replace in 2D Eq. (Bl) 
{2/y/TT) dse~^ ReJ with ^/n rfs erfc(s)Re J. As 
in 2D, the s independent part of ReJ enters the result 
being unchanged since -^7? erfc(s) = 1, whereas 
Vtt /o°° ds erfc(s) Ins = —1 — 7/2. This gives the 3D ver- 
sion of Eq. (23). 



it for the 2D isotropic case of Eq. (19). To this end, one 
separates the integration domain in two: < s < 
and < s < 00. In the first, the integration variable 

is changed to y = s/^Ju whereas in the second to y = 
{s^Ju)~^. Then we obtain: 



S{u, a) 



ITU Jo - (l + y^y 
exp I ) — cos(7rcr) exp 



(CI) 



APPENDIX C 



An explicitly real representation of the integral (16) is 
given in Ref. 10 for the Fermi sphere. Here we provide 



the form easy to deal with numerically. 
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